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Abstract 

We consider 1 + 1 dimensional SU(N) gauge theory coupled to a multiplet of massive Dirac 
fermions transforming in the adjoint representation of the gauge group. The only global symmetry 
of this theory is a U(l) associated with the conserved Dirac fermion number, and we study the 
theory at variable, non-zero densities. The high density limit is characterized by a deconfined Fermi 
surface state with Fermi wavevector equal to that of free gauge-charged fermions. Its low energy 
fluctuations are described by a coset conformal field theory with central charge c = (N 2 — l)/3 
and an emergent M = (2, 2) supersymmetry: the U (1) fermion number symmetry becomes an 
i?-symmetry. We determine the exact scaling dimensions of the operators associated with Friedel 
oscillations and pairing correlations. For N > 2, we find that the symmetries allow relevant 
perturbations to this state. We discuss aspects of the N — > oo limit, and its possible dual description 
in AdS% involving string theory or higher-spin gauge theory. We also discuss the low density limit 
of the theory by computing the low lying bound state spectrum of the large N gauge theory 
numerically at zero density, using discretized light cone quantization. 



* On leave from Department of Physics and Center for Theoretical Science, Princeton University. 
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I. INTRODUCTION 



An important aim of many applications of the AdS / CFT correspondence to condensed 
matter physics is the description of quantum matter at variable, non-zero densities. Here 
'density' refers to the conserved charge, Q, of a global U(l) symmetry of the underlying 
quantum field theory in d spatial dimensions. Our interest here will be restricted to zero 
temperature states which do not break translational symmetry or the global U (1) symmetry. 
Thus, we will not consider 'solids,' 'charge density waves' or 'superfluids.' In the traditional 
phases of condensed matter physics, the only remaining possibilities for non-zero density 
states are the Landau Fermi liquid in dimension d > 2, and the Luttinger liquid in dimension 
d — 1. Both these states are characterized by a Fermi momentum, hp, whose value obeys 
the Luttinger relation: the volume enclosed by the (d— 1) -dimensional surface in momentum 
space at hp is proportional to the density, Q/L d , with the same proportionality constant as 
that for free fermions (L is the spatial size which we will take to infinity). 

Any other realization of quantum matter whose density can be varied continuously by 
an applied chemical potential can generically be referred to as a 'strange metal.' A very 
promising candidate of a strange metal is a model of fermions at non-zero density coupled 
to an Abelian or non-Abelian gauge field. The non-Fermi liquid effects are strongest in 
d = 2, and this model has been the focus of much study in the condensed matter literature 
[THT3]. The theory scales to strong coupling, and a perturbative expansion in the gauge 
coupling constant cannot be used to analyze the leading infrared behavior. The flavor large 
Nf expansion also leads to difficulty: an expansion in the inverse number of fermion flavors 
cannot be reduced to counting fermion loops because of infrared divergences [9j [TO] . 

Another possible approach is to take the gauge-charged fermions in the adjoint represen- 
tation of the gauge group, and to then take the 't Hooft large N limit for the SU (N) gauge 
group. In this case, infrared divergences do not spoil the naive counting in powers of N, and 
so even the non-zero density case has a 1/N expansion controlled by the genus of the surface 
defined by a Feynman graph, as in all matrix models [H] . However, one is then left with the 
generally intractable task of summing all graphs with a given genus. For certain supersym- 
metric gauge theories, such matrix field theories can, in principle, be solved in the large N 
limit by the AdS/CFT correspondence [T5HT7] . Studies of such finite density models by the 
AdS/CFT correspondence have so far only provided a rather incomplete picture of 

the non-zero density quantum state. The boundary theory has density, Q/L d , which scales 
as N 2 , and essentially all of this density is associated in the bulk with degrees of freedom 
which are beyond the infrared horizon, with an unknown fate. Under appropriate parameter 
regimes, gauge-invariant probe fermions ('mesinos') can acquire a Fermi surface; however 
such a Fermi surface is only associated with a density of order unity, and is incidental to the 
physics of the non- Fermi liquid state [36H£T] . These probe Fermi surfaces are analogous to 
conduction electron Fermi surfaces in the 'fractionalized Fermi liquid' state of Kondo lattice 
[121 H3] , and do not yield much information on the underlying non- Fermi liquid state. In 
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certain uncontrolled computations, all of the boundary density Q/L d can be associated with 
visible Fermi surfaces in the bulk [HHI9], but then the resulting state is a Fermi liquid, 
although an interesting non-Fermi liquid state seems to have been obtained in recent work 

M- 

In an attempt to shed light on the difficult question of the fate of the 'hidden' matter of 
density proportional to iV 2 , this paper will examine the problem of adjoint Dirac fermions, 
at non-zero density, coupled to a SU(N) gauge field in 1 + 1 dimensions, i.e. for d — 1. We 
will show that a number of exact results can be obtained for general N, which we hope will 
help elucidate the structure of the large N limit of such matrix field theories. 

We consider the theory with Lagrangian 

C = Tr [# (i-fDp -m- ^7°) *] - -±-Tr F^F^ (1.1) 

z 9ym 

with a SU(N) gauge field A^, gauge field strength F^ V) gauge coupling gyM, and adjoint 
2-component complex Dirac fermions \I/ with mass m. The chemical potential /x couples to a 
global U(l) charge which is distinct from all the SU(N) gauge charges. Note that this U(l) 
is the only global symmetry of this Lagrangian. An analogous d — 1 model was examined 
earlier for adjoint Majorana fermions |50H5l] : in that case there is no global U(l) that can 
be coupled to a chemical potential; it was found that the ground state had an energy gap to 
all excitations, even at m = 0. As we will see here, just introducing a global U(l) by making 
the fermions complex is sufficient to transform the physics, and a gapless compressible state 
is obtained provided // is large enough, or when // = m — 0. 

This theory is characterized by three energy scales, m, QymVN, and fi. We will consider 
first the "high density" limit fi ^> m, gyuVN, where we can begin the analysis with a Fermi 
sea of free gauge-charged fermions. Next, we will consider the opposite "low density" limit, 
where m <C 9ym^/N while ji is comparable with m. Here we have to begin with an analysis 
of the spectrum of the SU(N) singlet excitations of the zero density vacuum: this will be 
carried out via the discrete light-cone quantization (DLCQ) [531 I57)] . 

We begin with a description of our results for the high density theory. The theory of free 
fermions has a Fermi wavevector related to the variable U(l) density Q/L by 

^ = (iV 2_i)^. (1.2) 

The iV 2 — 1 prefactor is a characteristic signature identifying this Fermi surface as that 
of gauge-charged fermions; this Fermi surface is 'hidden' [39] because the single fermion 
Green's function is not gauge-invariant. The Luttinger relation implies that this value of 
kp will not be renormalized [5Tj. We can analyze the infrared singular effects of the gauge 
interactions by writing down a continuum theory obtained by linearizing the fermions about 
the Fermi wavevector. Then, following a procedure standard in the condensed matter physics 
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FIG. 1: Energy dispersion of the Dirac fermions as a function of momentum. The full line is the 
zero of energy. The shaded region represents the occupied states. The filled circles are at ±k F . 



literature, we can express the Dirac fermion in terms of its right-moving and left-moving 
components at the Fermi surface 

*(x, t) ~ — i= (u(-k F )Mx, t)e tkFX + u(k F )ij L (x, t)e~ ikFX ) ; (1.3) 

V Itip 



where u(±k F ) are the standard Dirac basis spinors at the Fermi wavevectors; see (2.2) and 
Appendix [C] for a more detailed explanation, and see figure [T] for an illustration of this 
field redefinition. We will assume that ^l,r, and the gauge-fields, are slowly varying on 
the spatial scale kp 1 , and so all spatial integrals of fields multiplied by non-zero integral 
powers of e ±tkpX vanish. (In condensed matter models, such theories are obtained in the 
continuum limit of a lattice Hamiltonian, and in this context we are assuming that density 
is incommensurate, and so there is no 'umklapp' scattering.) An immediate consequence is 
that the resulting low energy theory has an emergent global U(l) conservation law: the total 
number of left-moving and right-moving fermions are separately conserved. We will denote 
these two £7(1) charges as Ql and Qr respectively, with Q = Q L + Q R . This U(l) x U(l) 
global symmetry will be crucial to our analysis. All operators appearing in the effective 
low-energy Lagrangian must have both Ql = and Qr = 0. 

As we will describe in section [XXJ, the high-density, low energy theory so obtained is a two 
dimensional conformal field theory (CFT), associated with the coset 

SU(N) N SU(N) N 

SU(N) 2N 1 • j 

of central charge 

N 2 -l 

c=^—. 1.5 



4 



The two dimensions of the CFT are the euclidean continuation of the original 1 + 1 di- 
mensions. With the requirement that the coset CFT have a global U(l) x U(l) symmetry, 
it actually has the M = (2,2) supersymmetry [5H]. For N = 2,3 the central charges are 
c = 1, 8/3, and then the theories coincide with the M = 2 superconformal minimal models 
[58] with c = 3k /(k + 2) for k — 1, 16 (in the k — 16 case we actually find a certain consistent 
truncation of the minimal model); the N > 4 theories were only briefly considered earlier 
[55] . For all N, the i?-charge symmetry of these M = (2,2) superconformal field theories 
(SCFTs) is U(l) x U(l), and this provides the needed global symmetry; the SCFT has no 
other global flavor symmetries. Note that this supersymmetry is an emergent symmetry at 
low energies and high densities; it is not a symmetry of the underlying Lagrangian. It is also 
remarkable that the diagonal .R-charge is conjugate to the chemical potential, /i, as has been 
assumed by fiat in many earlier higher dimensional studies of non-zero density quantum 
matter. 

These two dimensional SCFTs are our 'strange metals.' They are T = phases with 
variable density in models with only a global U(l) symmetry, but with a central charge 
which can become much greater than unity. The density fluctuations associated with the 
£7(1) symmetry cannot be represented by a gap less scalar field which is decoupled from all 
other sectors, as is the case for Luttinger liquids (exceptions for the N = 2, c = 1 case 
will be discussed in detail below). We note that the 'Bose metal' phases found in multi-leg 
ladder models in Refs. [60J also have c > 1 and only a global U(l) symmetry, although they 
are not expected to be described by our SCFTs. 

Armed with this construction of the SCFTs, we will compute exact scaling dimensions 
of gauge-invariant operators. An important observable which is sensitive to the presence of 
the underlying Fermi surface of the deconfined fermions is the Friedel oscillation in response 



to a localized perturbation coupling to the density. Upon perturbing the Lagrangian (1.1) 
via £ imp = C + A 5(x)p(x, t), where the density operator 

p = Tr(% *) = Tr U{^ L + ^ R + - e ~^^ R + ^^y^ , ( L6 ) 
V A 4 A 4 / 

the Friedel oscillation response is 

/ / w a cos(2k F x) 

(P(*)>imp « A |x 1 |2Af + • • • , (1.7) 

where is the scaling dimension of the operator Ty^i/j^i/jr), which we will call the Friedel 
operator in the CFT. Equivalently, we can relate the Friedel oscillation to an oscillatory 
term in the density-density correlator in the original system without an impurity: 

cos(2k F (x — x')) oN 
(p(x)p(x')} oc +■■■■ (1-8) 
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Our exact result for Ap, obtained by finding the smallest scaling dimension of operators 
with Ql = 1 and Qr = —1, is 



Ap = 1/3 for all N >2. 



(1.9) 



Observation of the oscillatory terms in (1.7) and (1.8) constitutes a measurement of the 



/cp in (1.2), and is a direct signature of the gauge-charged Fermi surfaces in these strange 



metals. Unfortunately, we do not determine the A" dependence of the missing proportionality 



constants in (1.7) and (1.8); the vanishing of this proportionality constant in the N — > oo 



limit is the presumed reason for the absence of such Friedel oscillations in existing studies 
via the AdS/CFT correspondence [61J. 

A second important observable is the fermion pair operator Tr (i/j^i/jr). Condensation 
of this operator leads to a superffuid ground state. The spatial decay of its two-point 
correlations is determined by its scaling dimension Ap. Our result for Ap was obtained by 
finding the smallest scaling dimension of operators with Ql = Qr — —1: 



A P = 1/3 for all N > 3. 



;i.io) 



For N = 2 there is no fermion pair operator Tt^l'Pr) m the CFT, and so in the original 
gauge theory we expect the two-point functions of Q = 2 and Q = — 2 operators to decay 
exponentially fast. Instead, the lowest CFT operator with Qi = Qr appears for N = 2 at 
Ql = Qr = —3 and has scaling dimension 3 (see section II A[ ). In Appendix |A| we review 
the Luttinger liquid of fermions with short-range interactions (e.g. the Thirring model at 
non-zero density), and find that it obeys Ap = 1/Ap; this identity is clearly not obeyed by 
the present adjoint matter theory. Appendix |A| also computes the values of Ap,p in models 
of fundamental Dirac fermions coupled to a SU(A^) gauge field. 



Finally, to assess the stability of the theory in (1.4) as a description of the low energy 



limit of (1.1), we have to determine the scaling dimensions of all perturbations allowed by 
symmetry: these are all operators with Ql = Qr = 0. Here we again find a distinction 
between N = 2 and AT > 3. For A^ > 3 the smallest scaling dimension of such an operator is 



dim 



Tr(*pl*p L *p R ip 



R) 



2(N -2) 
3iV ' 



(1.11) 



which is smaller than 2, and so always relevant. So the A^ > 3 SCFT2 is unstable to 
such a perturbation. We are not able to assess the A" dependence of the coefficient of such a 
perturbation, or the ultimate fate of the ground state. A natural conjecture is that this is an 
instability to a paired superfluid. In contrast, for N = 2 there are no relevant perturbations, 
and only a marginal perturbation. 

The low density limit will be considered in section |III[ Here we determine the spectrum 
of SU(N) singlet excitations above the zero density vacuum with the aim of using this as 
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input to describe the finite // state as a dilute gas of such states. We determine the mass 
M of the lightest state for a series of values of Q; we will obtain a dilute gas of such states 
for jj, > M/Q. So we need to determine the value of Q for which M/Q is a minimum. Our 
numerical analysis, carried out in the limit iV = oo, suggests that M/Q may accumulate 
to a dense set of decreasing values as Q becomes larger (see figure [3J. This suggests that 
M/Q becomes degenerate at some value in the limit of large Q (however, more extensive 
numerical work is needed to decide if the degeneracy is actually present). This degeneracy 
would indicate that even in the low density limit it is not appropriate to use a description of 
a dilute gas of gauge-neutral particles, and suggests the possibility that the gauge-charged- 
Fermi-sea description of the high density limit applies even at low densities. 



II. HIGH DENSITY 



In this section we will analyze the regime \i 2> m, gyMV N. We will derive the SCFT 2 in 



(1.4), and then analyze its properties in the subsequent subsections. 



We begin by writing the Hamiltonian for the free Dirac fermion in (1.1) in terms of 



particle, p, and hole, h, creation and annihilation operators introduced via (C6): 



Hn 



dk r 
2ti 



Tr 



[\Jk 2 + m 2 — ii)p\k)p{k) — (\/k 2 + m 2 



(2.1) 



where k is spatial momentum. This defines kp by /i = \jk\ + m 2 . Now we introduce the 
left and right movers by 



ipR(k) = p(k F + k) , ip L (k) = p(-k F + k), 



(2.2) 



and then linearize about kp by approximating 



Hn 



dk 
— 'i 
2ir 



-v k Tr 



(2.3) 



where the velocity v = kp/ \Jk 2 F + m 2 . We will henceforth set v — 1. Carrying out the same 
mapping to low energy degrees of freedom in the presence of the gauge field, we obtain the 
effective Lagrangian 



cff 



Tr 



^ R (d T - d x )^ R + ^{(d T + d x )^ L + (A T - A x ) [^L i) L \ + (A T + A x ) [4, i> R ) 



Tr F 



(2.4) 



Note that this theory is of the same form as (1.1), but after setting jj, — m — 0. The CFT 



structure of this theory becomes clearer upon writing the complex Dirac fields in terms of a 
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pair of Majorana fields ipl R , and a = 1,2: 

*Pl,r= ^(<r + (2.5) 
and then the Lagrangian becomes 

c cS = ^Tr [r R {d T - d x )r R + r L (d T + d x )r L + (a t - A a )r L r L + {a t + A x )r R r R ] 

\ — Tr F 2 . (2.6) 



2g 



YM 



As is well-known, each adjoint Majorana fermion is equivalent to a SU (N) WZW model at 
level N [62-66J, each with central charge (A^ 2 — l)/2. We assume the gauge theory is in the 
strong coupling limit (gyM oo), and then the integral over the gauge field reduces to a 
constraint: the vanishing of the currents Jl = ipi^L anc ^ Jr = ^r^r- ^ is easily verified 
that these currents obey a SU(N) Kac-Moody algebra at level 2N [67], and central charge 



2(A^ 2 — l)/3. The standard coset construction [63J then leads to the CFT in (1.4). 

The fact that the low energy CFT2 of the gauged adjoint Dirac fermion system has 
M = (2, 2) supersymmetry was demonstrated explicitly in (12) of [59]. This fact can also be 
seen to follow easily from an extension of earlier arguments. It is useful to write the CFT in 
(L4j) in the compact notation (N, N; 2N) as a special case of the general diagonal coset model 



(k,£; k + t) of SU(N), which is SU(N) k <g> SU(N) t /SU(N) k+e . Section 3 of [63] considered 
the coset model (N, £] N + £) and established that it had J\f = (1, 1) supersymmetry. We 
can obtain the second pair of supercharges by applying the same argument to the coset 
(k,N;k + N), and so conclude that the coset (N, N;2N) has M = (2,2) supersymmetry. 
This argument also shows that the i?-charge symmetry rotates between the two SU(N)n 



components i.e. between the two a components of the Majorana fermions. From (2.5) we 
then see that the -R-charge symmetry is the global U(l) which is conjugate to the chemical 
potential. 

In the following subsections we will describe the structure of these theories, including 
their modular-invariant partition functions and operator scaling dimensions. 



A. N = 2 

Let us first discuss the simplest non-trivial CFT, corresponding to the SU (2) gauge theory 
coupled to an adjoint Dirac fermion. This c = 1 CFT may be described by the coset 

SU(2) 2 ® SU(2) 2 

SU(2) 4 • [Z - () 
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The primary fields of this coset theory are therefore labeled by three SU(2) spins ji, J2, j 
where j = \j\ — J2I, • • • , ji + J2- Their conformal weights are 

fcW„*;i) = *^ + *^-^+». (2.8) 

The M = 2 super conformal symmetry fixes which values of (ji,j2,j) appear in the spectrum. 
Here n is a non-negative integer determined in terms of It will be zero for the 

cases of interest below. 

An important quantity characterizing a CFT2 is its modular invariant partition function 
on a torus: 

Z(r,f) = ^ e 27rir(fc > -c/24) e -27rif(fej-c/24) ^ ^.9) 
3 

where the sum runs over the entire spectrum, and (hj,hj) are the holomorphic and anti- 
holomorphic conformal weights of the state j. Once the modular invariant is known, it is 
not hard to read off the spectrum of the theory. While there is a variety of possible modular 
invariants at c = 1, the M = 2 superconformal invariance turns out to fix Z(t, f) completely, 
up to an additive constant. 

The c = 1 CFT turns out to be the simplest member of the series of M = 2 superconformal 
minimal models [58j EH] with central charges c = 3k /(k + 2): it is its k = 1 member. The 
dimensions of the M = 2 superconformal primary fields, and their £/(l) charges are in general 
given by [58] [59] 

p 2 — 1 — (s — r) 2 Irl s — r r , 

4(* + 2) + ¥' '=2(IT2) + 4- < 2 - 10 > 

where 1 < p < k + 1, \s\ < p — 1, and p — s must be odd. In the NS sector r = 0, while 
in the R sector r = ±1. The NS sector operators with p = s + 1 have h = q; these are the 
special operators that form the chiral ring [68]. We should note that each primary field of 
the extended M = 2 algebra gives rise to various Virasoro primary fields obtained by acting 
with the supercharges and U{1)r current oscillators. 

The k — 1 theory has the following NS-sector M = 2 primaries: the identity operator, 
and the operators with h — g and U(1)r charge q — ±g. In the R sector the primary fields 
are (h — ^,q — ±^) and (h — |,g = ±|). The former are the R ground states with 
h = c/24. 

Now, let us recall the c = 1 CFT of a compact massless scalar field <fi of radius r, i.e. <fi is 
identified with + 27rr. The torus partition function of such a theory has the simple explicit 
form 

oo oo 

Z scalaT (r,f) = \ V (r)\- 2 E e^e-^ fc H . (2.11) 
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The spectrum of left and right momenta is 



Tl Tl 

k L = - - wr /2 , k R = -+wr/2, (2.12) 

where n and w are the integer momentum and winding numbers, respectively. The left and 
right dimensions of the Virasoro primary fields, 

exp(ik L (p L + ik R (f) R ) , (2.13) 

are h = k\/2, h = k R /2. In addition, for a generic radius r, this CFT has certain primary 
fields with ki = k R = which occur for h = h = n 2 where n is an integer. The simplest of 
such discrete primary fields is the exactly marginal operator d<j)d<j) which changes the radius 
r. 

It is well-known that the compact scalar theory at radius r = 2a/3 (in units where the 
self-dual radius is \/2) has J\f = 2 supersymmetry. This theory is the bosonization of the 
above Af = 2 minimal model which provides the correct modular invariant. Let us consider 
some of the simplest operators in the bosonized theory and translate them into the original 
adjoint fermion language. The marginal operator that changes the radius, d(j)d<j), corresponds 
to Jl(z)J r (z), where 

J L oc Tt(V4,Vl) , Jr oc Tr(il>y R ) (2.14) 

are the U(l) x U(l) currents. We identify the U{1) x U(l) charges in the k — 1 N — 2 
minimal model as qi = and qR = — ^jj=. It follows that the relation between n and w in 
the compact boson model and the integer charges Ql and Qr of the fermions in the gauge 
theory is 

n = Q L -Q R , w = -(Q L + Q R )/6 , (2.15) 

where Ql = Q<1l and Qr = 6qR. 

For n = ±l,w = we get h = h = ^. This corresponds to h(h, §; 1) in (2.8). These 



two spin zero operators are products of the R-sector J\f = 2 superconformal primaries with 

For n = ±2,w = we get h n= ± 2 = h n= ± 2 = | corresponding to h(l,0; 1) or h(0, 1; 1). 
These two spin zero operators are products of the NS-sector M = 2 superconformal primaries 
with dimension |. In the gauge theory these operators are Tr^^^) with charges Ql = 
1, Qr — — 1; and Tt^iPl^r) with charges Ql = — 1, Qr = 1- Their sum is simply the fermion 
mass term. The total dimension of these operators is Ap = h + h = |; this is the scaling 
exponent for decay of the Friedel oscillations. 

For (n, w) = ±(3,0) we find h = h = |. The holomorphic part of this operator is the 
(h = |,g = ±|)A/"=2 superconformal primary from the R sector (in the coset theory, it is 

MU;o) = |). 

We could also consider n = ±4, w = operators with /i n= ±4 = h±4 = |. These operators 
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are not superconformal primary fields, but they are Virasoro primary. We note that h±4 = 
h±2 + This suggests that the n = ±4 operators are obtained from n = ±2 by acting with 
a holomorphic and an anti-holomorphic supercharge. 

For n = 0, w = ±1 we get h = h = |. According to (2.15), these are the operators 
with Ql = Qr = ±3 responsible for superconducting correlations. Since e iy ^"^ L is the 
supercurrent, we identify the n = 0,w = ±1 operators with double-trace operators which 
are products of two super-currents, each having dimension 3/2. Their net dimension is 3. 
We do not find, therefore, separate fermion pair operators Tr and Tr(i()' L if)' R ) . This 
is a special feature of the N = 2 case; we will see that such distinct operators appear for 
N>3. 



B. N 



In appendix [B] we outline a general approach to obtaining the modular invariant partition 
sum for the (N, N; 2N) cosets and, thereby, the operator content of the theory. We first 



show how for N = 2 this approach reproduces the result given in (2.11) and then proceed 
to the case N = 3. 

The proper starting point for the N = 2 gauged fermion model is the £0(6)1 invariant 
partition sum, broken down to SO(3)i x SO(3)±, 

Z SO(6)! _ 1 SO(6)!|2 + |^SO(6)i |2 + 2 |^SO(6)i |2 

= |xf° (3)l xf lo(3)l + X V 50( ^X V 5 ° (3)1 | 2 + 4|xf 0(3)l X v 50(3)l | 2 + 2|xfp 0(3)l xfp° (3h | 2 • (2.16) 

We consider the branching rules of the relevant combinations of SO(3)i characters. The 
M = (2, 2) superconformal symmetry of the coset CFT2 guarantees that chiral branching 
functions will organize into characters of M = 2 SCFT, defined as 

ch R,NS = ^^(Lo-c/24)] (2 17) 

(the complete characters of the M = 2 superconformal symmetry would keep track the U(l) 
charges as well - for simplicity we suppress those in our notations). Considering the explicit 
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leading terms in the various characters, we established the following relations 



50(3)1 50(3)1 



* -xi ■ >1 +x s v 0{3h x s v 0{3h + 2xT-x;' 



SO(3)i SO(3h 



CJ1 U(0) + X(4) J + Z cn l/6 X(2) 



SC/(2) 4 



NS SU(2) 4 



50(3)i 50(3)i 50(3)i 50(3)i 
X.1 X.1 "r a.v A-v 



ch o S (xf ^ (2)4 



50(3)i 50(3)i 
Asp Asp 

,R 



cn i/24 U(0) + ^(4) ) + cn 3/8 X( 2 



517(2)4 « 
A (4) 



5(7(2) 4 



2ch 



2xf O(3)i xf O(3)i 
NS SU(2) 4 



1/6 X( 2 ) 



5t/(2) 4 



(2.18) 



~ NS 

Note that the characters ch are obtained by inserting a (— 1) F , with F the fermion parity 

~ R 

operator. In the R-sector the characters ch are constants known as the Witten index of 
the sector. 

For SU{2) i there exists an exceptional invariant (labeled as -D4 in the literature |69j), 



which groups the characters according to the automorphism (B7) 



Z SU(2) 4 



|X( 



0) 



SC/(2) 4 |2 
+ X(4) 



2|X( 



51/(4)4,2 

2) 



(2.19) 



This invariant sets a modular invariant metric on the SU{2) i characters. Using the metric 
we can project out the level-4 characters in the various product terms in the SO(6)\ partition 
sum. The final result is the following modular invariant coset partition sum 



r7-coset(2,2;4) 



[|chf | 2 + 2|ch^ s 6 | 2 + (ch NS ^ch^)] 
i[|ch*| 2 + 2|ch R '-'i 



1/241 



(2.20) 



Comparing with the expression (2.11), evaluated at the J\f = 2 radius r = 2y / 3, one checks 



that the two expressions agree up to a constant which we can write as the sum of the R-sector 
Witten indices and which is by itself modular invariant 



2-coset(2,2;4) _ g 



1 



scalar 



2|ch x 



/24l 



scalar 



1 . 



(2.21) 



We are now ready to take on the case N = 3. The central charge c = 8/3 corresponds 
to the k — 16 entry in the minimal series of unitary M = (2,2) SCFT2. One thus expects 
that the partition sum can be expressed as a finite sum of terms of the form ch^~ 2 ch h / g / . 
Modular invariant partition sums of this type have been completely classified [70] - our 
challenge is thus to identify the correct entry from the (rather extensive) list. 
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The N = 3 strange metal starts from 16 fermions with partition sum 



750(16)i = I SO(16)i|2 1 I SO(16)i |2 1 01 SO(16)i |2 
IXl I ' lAv I ' IXsp I 

= Ixf° (8)l xf° (8)l + X V 50(8)1 X V 5 ° (8)1 | 2 + 4|xf° (8)l X v 50(8)l | 2 + 8|x s 5 p 0(8)l Xs 5 p 0(8)l | 2 • (2-22) 

A curiosity specific to iV = 3 is that the vector representation of 5*0(8) is isomorphic to the 
spinors - in the coset theory this leads to degeneracies between NS and R sector dimensions. 

As before we now study the branching into characters of SU(3)q times branching functions 
which we expect to be characters of a M^-algebra extension of M = (2, 2) superconformal 
symmetry. 1 Working through explicit details, one observes that in the r.h.s. of the branching 
rules, the SU(3)q characters consistently show up in the combinations 

SU(3) 6 , SU(3) 6 , SU(3) 6 SU(3) 6 , SU(3) 6 SU(3) 6 
X(qo) + A(60) * + " ^(06) ' *(11) + X(41) + X(u) ' 

51/(3)6 , SU(3) 6 SU{3) 6 SU(3) 6 (t) 9 ^ 

A(33) A(30) A(03) ' A(22) ■ ^.ZO; 

A modular invariant projector for these terms is provided by the Vq invariant of SU(3)q [H] 

Z SU(3) a _ | X(oo)+X(6o)+X(o6) |2 + | X(n)+X(4i)+X(i4) |2 + | X(33)+X(3o)+X(o3) |^ ^.24) 

Completing this analysis, we have identified (up to a constant) the partition sum for the 
(3, 3; 6) strange metal with the exceptional invariant of M = (2, 2) superconformal symmetry 
at c = 8/3 which in the classification |70j is labeled as M 4,2 with parameters v = 3, z = 1, 
x = 1. This invariant involves a subset of all NS and R characters at k = 16; furthermore, 
those that survive are grouped into groups of 6 (4 groups in each sector) and per sector 
one group of 3. So in total, 54 fields survive. In the NS sector the extended characters are 
(A/" = 2 characters labeled as ch; =p _i iS with r = for NS and r = — 1 for R) 

h NS,ext _ ,NS , h NS , iNS , jNS , jNS , ijVS 
cn — ui -1- cn 16 -|- cn 166 -|- ui 16 _ 6 ui 1612 -|- ui 16 _ 12 

h NS,ext _ ,NS , h NS , h NS , jNS , jNS , jNS 
U1 l/9 ~~ LI1 2,0 ^ Ln 14,0 ^ L11 14,6 1^ L11 14,-6 ^ L11 14,12 ^ ul 14,-12 

h NS,ext _ ,NS , h NS , jNS , ,NS , jNS , jNS 
U1 l/3 ~~ U1 4,0 ^ LI1 12,0 ^ U1 12,6 ^ ul 12,-6 ^ ul 12,12 ^ ul 12.-12 

= < + <o + <*$ + chf„ 6 + < 6 + <_ 6 
</Y? = ch™ + < + chf_ 6 (2.25) 



1 See [71] for a general review of W-symmetry and [72l [73] for VF-extensions of J\f = 1 superconformal 
symmetry in general cosets involving a SU(N) n factor. 
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and in the R sector 



i K,ext _ , K , . K , i K , i K , An , i K 

U1 1 _ U1 16,-16 + U1 16,14 + U1 16,-10 + Ln 16,8 + ul 16 -4 + ul 16,2 

i R.ext _ i R , i R , i R , i R , i R , i R 

1/9 ~~ Lli 14,14 ' L,Ii 14,-10 ' ul 14,8 ' L ' Ii 14,— 4 ^ LIi 14,2 ^ L11 2,2 

i Rext _ i R , i R , , R , i R , , R , , R 

1/3 ~~ LIi 12,-10 L,11 12,8 L,Ii 12,-4 UJ 12,2 ^ LIi 4,-4 ' L11 4,2 

i R,ext _ , R i R , R i R , R , R 

U1 2/3 ~~ U1 10,-10 ^ U1 10,8 ^ LI1 10,-4 ^ ul 10,2 ^ LI1 6,-4 ^ LI1 6,2 

ch^f = ch*_ 4 + ch* 2 + ch* 8 . (2.26) 

The extended vacuum character has the following content (returning to the notation ch^ j? 
for the M = 2 characters) 



ch, 



NS,ext 




(2.27) 



L11 /i=0,g=0 ^ L11 /i=2,g=l/3 ^ L11 fe=2,g=-l/3 



+ cn /i=7/2, g =l/6 + Cn fc 



NS i r l,NS 

7/2,<j=-l/6 Lli / l =4,g=0 



These fields span a P^-algebra extension of M = 2 superconformal symmetry, with extra 
currents of dimension 2, 2, 7/2, 7/2, 4. The IV-currents at h — 2 are given by 



W L = T^M) , ^2 = Tr(^f) 



(2.28) 



These currents exist for all N > 3. For N = 2 the triple products ^| and ip^ are SU{2) 
singlets and the trace vanishes. For N > 3 explicit expressions for Wl and W\ in terms of 
component fields {iPl)a an d {Wl)a-> with A = 1, 2, . . . iV 2 — 1 an adjoint index, involve the 
3-index (i-symbols cIabc- 
The M 4,2 invariant reads 



E 



Ich 



x ^ x f + 2|chf 1 s ; 1 e fi 2 +(ch NS ^ch 



NS. 



h=0,l/9, 1/6,1/3 



+ 



E 



Ich 



R,ext 1 2 
ft 



R,ext 1 2 



+ 2|chf/9< 



and the claim is 



h=l/9, 1/3,2/3,1 



^coset(3,3;6) rr 



M 4 > ; 



1/3. 



(2.29) 



(2.30) 



The U(1)r charges of the fields surviving in the partition sum are easily extracted from the 
field labels in the partition sum. In both the NS and the R sectors, the U(1)r charges are 
multiples of ±1/6. 

Comparing with N = 2, we observe, at iV = 3, the presence of fermion pair operators 
Tr^ipLipji) and Tr^j^Jj) which were absent from the operator spectrum for N = 2. From 



(2.25) we read off that these operators are in the same extended symmetry multiplet as the 
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current 


{Ql, Qr) 


(M) 




J L oc Tr(^ L ^I) 


(0,0) 


(1,0) 


N > 2 


Gl oc J-r(^2J 


( 1 


(2,0) 


iv > z 


G L oc Tr(C) 


(§,o) 


(|,o) 


iV > 2 


T L oc Tr^L&ili + ili&if> L ) 


(0,0) 


(2,0) 


N > 2 


W L ex Tr(V[V!) 


(-i.o) 


(2,0) 


N > 3 


W L ex Tr(i/> L 1$) 


(5.0) 


(2,0) 


N > 3 



TABLE I: This table lists some of the (left) chiral fields (currents) in the (N, N; 2N) coset model. 
The currents Jl,Gl,G^ l ,Tl constitute an M = 2 superconformal algebra. The dimension-2 cur- 
rents Wl, Wt are the first of an extensive set of extra primary symmetry generators that exist for 
N > 3. 



mass (Friedel) operators Tt(^ l ^r) and Tr^i/^Jj). Both these sets of operators have scaling 
dimension 1/3 - the degeneracy is due to the fact that the zero mode of the charged spin 
two currents in ( |2.28 ) relates the two sets of operators. 

We also observe the presence of a number of charge neutral relevant operators. The most 
relevant operator is the (I — 2, s — 0) field in the NS sector, with scaling dimension A = 2/9, 
as in (1.11). 



C. General N 

We have also obtained several explicit results for N = 4 and higher, confirming the 
general structure outlined above and in appendix [Bj We defer the details of the description 
of the general A" case to a forthcoming publication. Here we briefly indicate some of the 
findings, paying special attention to the extrapolation to large N. 

• In tables [T] and [IT] we list the general form of the leading currents and primary field 
operators in the (N, N; 2N) coset. These operators are all in the NS sector. We see that 
their scaling dimensions are either ^-independent or are such that h (h) converge in 
the large A" limit to nuR\/Q if the operator involves tiur\ left (right) moving fermions. 
In particular, we see that both fermion pair operators and the Friedel operators are 
present in the spectrum and have a scaling dimension Ap = Ap = 1/3 independent 



of N. Their degeneracy is due to the presence of the charged spin two currents (2.28) 
which are part of the extended jy-symmetry for any N > 3. 
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operator type 




(h, h) 


channel 






i\ i\ 

V6' &> 




(10. 


.01) 


JM > 2 


Tr(^) 


V 6' 6^ 


V6' 6^ 


(10. 


.01) 


N>2 


Tr(^L^i?) 


/ 1 1 \ 

V 6' _ 6^ 


/I 1 \ 

(I'D 


/in 

(10. 


.01) 


N > 3 


Tr(^) 


/ 1 1 \ 

U' 6 J 


(I'D 


/in 

(10. 


.01) 


N >3 




r 1 M 

V 3' 3^ 


U> 3^ 


(20 . 


10) 


AT ~> 1 
iv ^_ 


Tr(V'iV'iV'iiV'fl) 




(3' 3) 


(20 . 


10) 


AT ^> Q 


TV./,/, „/, „/,t „/.t \ 


( 1 M 

V 3' 3^ 


U> 3^ 


(20 . 


10) 


/\7" "~> Q 


TrfolJ nil nil nil \ 




(3' 3) 


(20 . 


10) 


iv ^ 


TV/V,! „/■ „/,t „/, N 
Lt(1p L 1p L 1p R 1p R ) 


(U,U) 


U> 3^ 


(20 . 


10) 


iv 




(0, 3) 


U> 3) 


/on 


1 r\\ 
• iU J 


AT ^ Q 
iv ^ 




/i 1 \ 

V3' ~~ 3/ 


/ 1 In 

(I'D 


/on 
{20 . 


■ 10) 


iV > 3 


m / / T / T / T / \ 


(1°) 


(DD 


(20 . 


.10) 


> 3 


4-4-4-4- 

Tr(VM^) 


(- ~) 
V3' 3^ 


(DD 


(20 . 


.10) 


N > 3 


same 9 terms 




(DD 


(010. 


..02) 


> 3 


same 9 terms 




(- -) 

V3' 3^ 


(10. 


.01) 


AT > 3 




(0,0) 


/JV-2 JV-2\ 
v 3N ' 3/V / 


(20. 


.02) 


> 3 




(0,0) 


/ iV+2 JV+2 \ 
v 3N ' 3/V / 


(010. 


.010) 


A^ > 4 




(_1 1) 

V 6' 


/JV-2 JV-2\ 
V 2/V ' 2N > 


(30.. 


.011) 


A^ > 4 



Tr^lMy R ) 2 (0,0) (^,^P) (210... 012) iV > 4 
Tr{1>lMy R ) 2 (0,0) (^,^) (40... 020) iV > 4 



Tr(^ L ^ R ) n (0, 0) (n^p, n^=^) (2n0 . . . 02n) N > n + 2 



TABLE II: Operators in the (N, N; 2N) coset model. The operators listed (all in the NS sector) 
are primaries with respect to the N = (2, 2) superconformal algebra. The channel indicated is the 
SU(N) representation of the gauged diagonal SU(N)2N subalgebra. We give the complete list of 
2- and 4-fermion primaries and in addition list some of the higher operators. 

• One can also quickly verify that the scaling dimensions in the R sector satisfy 

R c A^ 2 — 1 , 
h R > — = —— 2.31 
~ 24 72 y 1 
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so that for large N all R-sector operators are irrelevant. 

The algebraic structure of the general-A" coset appears to be highly involved. The 
cosets of the form 

SU(N) k+1 {2 - 62) 

are known to have an extended Wn symmetry with one chiral current for each spin 
s = 2, . . . N (TTJ [75] . For a general diagonal coset (which includes our coset as a special 
case) 

SU(N) k ® SU(N) e 

SU(N) k+e {2 - 6S) 

we can construct higher spin-s currents from polynomial (of degree s) combinations of 
the individual numerator SU(N) currents which commute with the diagonal SU(N). 
See (7.42) and (7.43) of [71] for an explicit expression for the spin-3 current. These 
general cosets are expected to have many additional currents as well. In the case of 



N = 3 we found two charged spin two fields, see (2.28), in addition to the stress tensor. 
Since these are primaries of the M = 2 supersymmetry, this then implies the existence 
of charged spin | and spin-3 currents as well. Similarly, we see from (2.27) that there 
are charged spin | currents and their Af = 2 descendants. 



For general N we expect to be able to build many chiral operators from ip^, i\r L and 
holomorphic derivatives. In the large N limit we do not expect any trace relations for 
such operators of spin s <C N. Then by the usual counting arguments for words built 
from matrix valued fields one expects to roughly see a Hagedorn growth in the number 
of these currents (as a function of the dimension, which is the same as the spin). Thus 
we might expect an algebra much larger than the conventional Wn symmetry algebra. 
It would be interesting to work out the consequences of this larger symmetry algebra. 

There are a number of issues that would be important to explore further for this class of 
coset models. 

A key algebraic structure of M = 2 SCFT2 is the chiral ring: the collection of NS sector 
primaries with q = ±h, which form a closed algebra under fusion [6H]- In the N = 3 case, 



we read off from the NS sector primaries in (2.25) that the only chiral primaries are those 



with p — 1 = s — 6, 12 in (2.10). The chiral ring in this minimal model case is generated by 
one generator x obeying the relation x k+2 = x 18 = 0. The chiral primaries present in the 
coset correspond to the elements x 6 and x 12 which form a consistent subring of the original 
chiral ring. The large N structure of the chiral ring can, we believe, be exploited to study 
the large N physical characteristics of the gauged fermion model. 

Finally, it would be very interesting to uncover the AdS dual to this interesting class of 



CFTs (in the large N limit). The simpler coset models in (2.32) have been identified [76l[77j 
with a class of higher spin Vasiliev theories on AdS% [7S]. These have a single bulk gauge 
field of every spin s > 2. Roughly speaking there is a single Regge trajectory in the bulk 
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and no Hagedorn growth of states. This is related to the fact, mentioned above, that there 
is a single conserved current in these coset CFTs for a given spin. The jy-symmetry in 
the coset considered here is very different having many more fields. In general, we expect 
there to be a Hagedorn density of states in this theory corresponding to single trace words 
built from ipL^^ipR,^ together with derivatives (modulo the projection by the diagonal 
SU(N) currents). Thus the bulk dual is presumably a string theory on AdS^. This fits with 
the fact that these cosets have a central charge proportional to N 2 as opposed to the cosets 
in (2.32) which have c oc N and behave like vector models with fewer gauge invariant states 
[79] [76]. 

Though the coset theory is a strongly interacting CFT as evidenced by the anomalous 
dimensions of different operators, it is interesting that the dimensions of the operators in 
the table [IT] are like those of a free theory in the large N limit. They are integer multiples 
of |. In that sense the spectrum is similar to that of free Yang-Mills theory or that of the 
D1-D5 CFT [80l El] at a symmetric orbifold point. Note that these are the points where 
we expect a dual tensionless string theory with an unbroken higher spin symmetry. This is 
consistent with the fact that these cosets have a large chiral iy-algebra as discussed above. 



III. LOW DENSITY 

In this section, we will study the behavior of the SU(N) gauge theory coupled to two 
adjoint multiplets of Majorana fermions for small U(l) chemical potential /i. We will work 
in the regime where m, the mass of the adjoint, is much smaller than the scale set by the 
't Hooft coupling gyM^/N- We are interested in the low energy dynamics of the system as 
H is increased from zero to values comparable to the scale set by m. Some features of this 
low energy theory can be inferred from the spectrum of gauge invariant hadronic states at 
(jl = 0. Let us describe how this works for our model. 

Generally, it is extremely difficult to compute the spectrum of hadronic bound states in 
gauge field theories such as QCD. In 1 + 1 dimensions, however, the light cone quantiza- 
tion can make this problem tractable. Compactification of the light-like coordinate on a 
circle: x~ ~ x~ + L, a formal regularization procedure known as the discrete light cone 
quantization (DLCQ) [55j|56], typically reduces the problem to matrix diagonalization. The 
problem is further simplified in the planar (large N) limit where the SU(N) singlet states are 
non-interacting. The physical bound state spectrum can be inferred by taking the decom- 
pactification limit for the light-cone coordinate x~ . This can be computationally expensive 
for certain models, but is nonetheless a well controlled approximation scheme. Often, this 
continuum limit is presented by taking the limit K — > oo, where the integer 'harmonic 
resolution parameter' K enters the definition light cone momentum: 

. , iiK 

P + = — (3.1) 
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After diagonalization of the light cone hamiltonian, P , the spectrum of bound state masses 
M is read off from 

M 2 = 2P + P~ . (3.2) 

In this paper we will be primarily concerned with the spectrum of bound state masses 
in the regime where m 2 is kept fixed while the 't Hooft coupling gy M N is sent to infinity. 
We find a rich spectrum of bound states whose masses divided by m approach constants in 
this limit. These bound states originate from the large iV coset CFT discussed in section II, 
perturbed by the operator m(Tr(^-^] ? ) + h.c.) of dimension 1/3. This mass operator breaks 
the Ql — Qr symmetry of the CFT, but the overall U(l) charge symmetry, Q = Q L + Q R , 
remains unbroken. Therefore, the bound states break up into sectors labeled by the integer 
charge Q. 

The DLCQ of a closely related system consisting of a single adjoint multiplet of Majorana 
fermions has been analyzed in the past [5"U| 152"! |5"1] . However, in the same limit of sending 
the 't Hooft coupling to infinity while keeping m 2 fixed, all bound states acquire masses of 
order ~ gYM^/N . The fact that there are no light bound states with masses of order m is 
due to the triviality of the CFT arising from the SU (N) gauge theory coupled to an adjoint 



Majorana fermion: in that case, instead of (1.5) one finds c = 



We will summarize the details of the DLCQ computation in Appendix C. In this section, 
we will focus on presenting the results of the computation. 

First, let us describe the basic physics of the model. In 1 + 1 dimensions, gauge fields are 
non-dynamical but serve as agents binding the colored matter fields. The hadronic bound 
states will then be superpositions of traces of products of the adjoint creation operators 



(C22). The problem can be separated into boson and fermion sectors depending on whether 
the number of the creation operators in the trace is odd or even. 2 

When there are two adjoint Majorana fermions, there is a U(l) = SO(2) global flavor 
symmetry which provides an additional quantum number to label the states in the spectrum. 
One way to make this manifest is to combine the two adjoint Majorana fermions into an 
adjoint Dirac fermion, and count the difference between the number of fermions and anti- 
fermions for each state. 

Suppose for /i = we succeed in computing the masses M and charges Q of all the bound 
states. Suppose also that all states have M > and as a result the theory is gapped in the 
far IR. What happens when /i is increased? 



2 The CFT arising in the m = limit exhibits the J\f = 2 supersymmetry relating the bosonic and the 
fermionic operators. Since the light cone quantization is unreliable in the presence of massless states, we 
will keep m non- vanishing. Then the supersymmetry is broken; so, the boson and fermion bound state 
spectra are not identical. However, the bound states with m <^ M gym^N may exhibit approximate 
supersymmetry. Such highly excited states are difficult to access numerically, but it would be very 
interesting to look for the emergent supersymmetry in this region of bound state masses. 
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We expect that some of the particles will condense when 

( = M-fj,Q (3.3) 

becomes negative. 

If we were working in dimensions greater than 1 + 1, then depending on whether the first 
state for which ( becomes negative is a boson or a fermion, the system would exhibit the 
universal behavior of either the Bose-Einstein condensation or formation of a degenerate 
Fermi gas. In 1 + 1 dimensions, the distinction is somewhat blurred by the fact that bosons 
and fermions are related by bosonization. If the spectrum is sufficiently generic so that there 
is precisely one state for which ( is going to zero at the minimal critical /i, then one expects 
the system to behave as a Luttinger liquid. 

Additional subtleties can arise from the fact that there might be a degeneracy causing 
more than one state to hit ( = at the same time. Logically, there are three distinct 
possibilities. 

I The value of ( goes to zero for exactly one state. 
II The value of ( goes to zero for several, but a finite number of states. 
Ill The value of ( goes to zero for an infinite set of degenerate states. 

The DLCQ computation of the bound state spectrum can help distinguish among these 
three possibilities. 

The details regarding the implementation of the DLCQ procedure are summarized in the 
appendix. Here we will merely present the result, where we display the full spectrum in 
figure [2] 

The points illustrated in figure [2] are to be interpreted as follows. 

1. These points correspond to the spectrum of fermions for which the harmonic resolution 
parameter K takes odd integer values. One expects to recover the exact spectrum in 
the limit K — > oo. 

2. The 't Hooft coupling gy M iV is taken, for the sake of definitiveness, to be 2tt ■ 10 3 
times the bare mass-squared of the fermions, m 2 . As long as this number is very large, 
the spectrum in the range illustrated, presented in units where m? = 1, is insensitive 
to its precise value. The idea is to extract the behavior of this system in the limit of 
large gy M N. 

3. M is the mass of the hadronic bound state. Here we are displaying M/Q. The state 
with lowest Mj Q is the one whose ( will hit zero first as /i is increased. 

4. For each K, Q in the range from Q — 1 to Q — K are possible. However, the Q = 1 
states are heavier than the range of M illustrated in figure [2] and as a result are not 
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FIG. 2: The spectrum of fermionic bound states for K = 5, 7, 9, 11, 13, 15, 17, 19. The states with 
the same K are shown using the same color. Increasing K for fixed Q is illustrated by a gradual 
shift to the right in each of the columns. 

displayed in this figure. The state in the Q = K sector is decoupled from the rest of the 
dynamics. In general, taking the large K limit with Q fixed will give rise to a reliable 
extrapolation of the spectrum of that fixed Q sector. The spectrum with Q of the order 
of K, however, is sensitive to the DLCQ artifacts and does not effectively approximate 
the spectrum in the continuum limit. Here we have computed and presented the states 
for Q in the range 3 < Q < K — 2 with the exception of some small Q states for large 
values of K for which the computations became numerically intractable. 

5. For each value of K, states with different charges Q are displayed in separate columns. 
States with the same K, however, can be identified by the fact that they are plotted 
using the same color. For each Q, increase in K is indicated by gradual shift in the 
column of points to the right. 

In order to identify the states with smallest M/Q, one must, for each Q, track the lowest 
mass state and extrapolate to large values of K. The spectrum illustrated in figure [2] suggests 
that, for each Q, the masses are gradually increasing in a similar manner as K is increased. 
It also suggests that these increasing masses are converging to the large K limit. 

In order to analyze the asymptotic behavior of the masses of the low-lying states, it is 
useful to plot their masses as a function of 1/K. Since we are interested in how these states 
are affected by the chemical potential, we will actually plot M/Q as a function of 1/K for 
different values of K. This is illustrated in figure [3] for Q = 3, 5, 7, 9, 11. 

We have also superposed a line indicating the linear extrapolation available from the set 
of data available. These lines cross 1/K = at finite values. This can be viewed as a crude 
method to extract the extrapolated value for the large K limit. 
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FIG. 3: The spectrum M/Q of lightest states as a function of 1/K for Q = 3,5,7,9,11. The 
lines are linear fits to the available data. The points for Q = 3 have the largest M/Q. For each 
successive Q, the M/Q is decreasing. The intercept of the linear fit at 1/K = is the extrapolated 
value of the M/Q for the lightest state for fixed value of Q. 

What we see in figure [3] is the tendency for the large K limit of the Mj Q of the lightest 
states to become dense with increasing Q. Since Q can get arbitrarily large, this suggests 
that in the large K limit, states with arbitrarily large Q are converging to the same value of 
M/Q. In other words, we seem to be finding out that our model is exhibiting scenario III 
enumerated earlier in this section. It should be kept in mind, however, that one must send 
K — > oo first, and then look for a trend as Q is increased. The order of these two limits 
cannot be exchanged. 

It is not easy to determine the critical value fi cr u of M/Q. To settle this question, a higher 
precision computation is required. Unfortunately, for the reason given above, one must 
explore very large values of K in order to explore large values of Q. This is computationally 
very expensive. 

The fact that there may be infinitely many states with degenerate fi cr i t = M/ Q suggests 
that, when \i approaches fi cr u, the system may undergo a transition from a gapped phase 
into a phase with a non-trivial interacting conformal field theory with c > 1 in the IR. 

The DLCQ results presented in this section focused mainly on fermionic bound states. 
We have also carried out the analogous computation for bosonic bound states and found a 
similar behavior of M/Q. We will not present the detailed results of our analysis for the 
bosonic spectrum here. Some technical issues which arise for the DLCQ computation in the 
bosonic sector will be described briefly in Appendix C. 

Let us, in closing, mention that since all of these results were presented for the case 
where N was taken to infinity first, it is certainly possible for the seemingly dense spectrum 



22 



of M/Q to be discretized by a fine structure of order N~ a for some a > 0. Such a structure 
can change the basic feature of our model from scenario III to II or I. This may also be 
related to the fate, following the relevant flow, e.g. (1.11 ), of the coset CFT discussed in the 
earlier sections. 



IV. CONCLUSIONS 

In this paper, we described a novel metallic state of matter in one spatial dimension, with 
a continuously variable density. 

The well-studied one-dimensional metallic state is the Luttinger liquid, and in many 
respects this state can be considered to be the natural limiting case of the Fermi liquid 
state of higher dimensions. Indeed, the Luttinger liquid reduces to a free fermion model 
at a specific value of the Luttinger parameter, and other parameter values are continuously 
connected to this one. It has central charge c = 1, and this is directly linked to the massless 
scalar representing fluctuations of the globally conserved U(l) charge density. We note 
that some Luttinger liquids have additional gapless modes associated with other global 
symmetries (and so a larger central charge), such as the models described in Appendix [A} 

Our state was obtained by considering a non-zero density of Dirac fermions carrying 
adjoint color charges of a SU(N) gauge field. We found a 'strange metal' state described 
by two-dimensional superconformal field theories with central charges c = (A^ 2 — l)/3. For 
large A" the central charge becomes large, while the global symmetry remains only the U(l) 
associated with fermion number conservation. The strange metal has a Fermi surface with 



a Fermi wavevector, given by (1.2), which is equal to that of non-interacting color-charged 
particles; the Fermi wavevector changes continuously as the density is varied. This Fermi 
surface is 'hidden' , because the single fermion Green's function is not a gauge-invariant 
observable. Nevertheless, the Fermi surface and the value of the Fermi wavevector are 



detectable in the Friedel oscillations of (1.7). We propose that this variable density state, 
with its large phase space of low energy excitations linked to its large central charge and 
Fermi surface of color-charged fermions, can serve as a paradigm for non-Fermi liquid states 
in two and higher spatial dimensions. 

The structure of our d — 1 strange metal is quite analogous to the 'hidden' Fermi surface 
states obtained recently for general d in Refs. SI] via the AdS/CFT correspondence. 

In Refs. [38, [39] it was postulated that the Fermi surfaces of gauge-charged particles could be 
detected by the hyperscaling violation of the thermal entropy density, and by a logarithmic 
violation of the 'boundary law' of entanglement entropy. The coefficient of the entanglement 
entropy logarithm was used to deduce a value for the Fermi wavevector which depended upon 
ultraviolet details only through the value of the density Q/L d , in just the manner expected 
from the Luttinger relation |39j. In our d — 1 strange metal here, analogous properties of the 
entropy and entanglement entropy are trivially satisfied, because the hyperscaling violation 
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exponent 9 = d — 1 = (in the notation of Ref. [39]), and every CFT2 has a logarithmic 
violation of the boundary law of entanglement entropy j82j [83]. However, here we were 
able to detect the Fermi surface, and determine the Fermi wavevector, from the Friedel 
oscillations in the density correlations. Obtaining the Friedel oscillations in the AdS/CFT 
description of strange metals in general d is clearly one of the important challenges for 
future work. In this direction, it would be useful to understand the large N dependence of 
the proportionality constant in (1.7): this should shed light on how Friedel oscillations of 
gauge-charged particles appear in holographic theories. We also note the interesting recent 
computation of [H] which detected Friedel oscillations in an anisotropic quantum liquid of 
strings in 5 + 1 dimensions. 

In section III we studied the bound state spectrum of the large N gauge theory in the 
limit 9ymN 3> m 2 using the numerical DLCQ approximation. This approach sheds light on 
the properties of the SU(N) gauge theory at low density. A very interesting phenomenon 
that we have uncovered is the emergent Af = 2 supersymmetry of the gauge theory in 
the limit m — > 0. A useful direction for future work would be to obtain some numerical 
DLCQ evidence for the emergent supersymmetry by studying the masses of highly excited 
bound states with gyM^/N 3> M 3> m. Hopefully, this spectrum will exhibit approximate 
supersymmetry. 

Another intriguing direction for future work, which was discussed at the end of section 
II. C, is the possibility of a dual description of the large N CFT in terms of a theory with 
higher-spin gauge symmetry in AdS^. The existence of such a dual description is suggested 
by the large ly-symmetry, and by the fact all the operator dimensions appear to approach 
constants in the large N limit that are quantized in units of 1/3. A useful 3 + 1 dimensional 
analogue of the Af = 2 supersymmetric large iV CFT we are considering may be the M = 4 
supersymmetric SU (N) gauge theory, which is dual to the AdS§ x S 5 background of type IIB 
string theory [T5TH7] . Both theories have anomaly coefficients of order N 2 . When the Af = 4 
gauge theory has a large 't Hooft coupling g YM N, the dual AdS 5 x S 5 background becomes 
weakly curved. In this limit the dimensions of all operators not protected by supersymmetry 
become very large. On the other hand, at vanishing 't Hooft coupling all the operators have 
integer dimensions, and their number exhibits the exponential Hagedorn growth. In this 
limit the curvature of the dual string background becomes very large; this is often referred 
to as the "tensionless string limit." It has been suggested [M] that a useful dual description 
of the free Af = 4 SYM theory may involve higher-spin gauge theory in AdS§ [85] coupled 
to an infinite number of additional fields. 

Similarly, if an AdS 3 string dual of our Af = 2 supersymmetric coset CFT is found, 
we expect it to be strongly curved. It is therefore interesting to ask if the coset CFT has 
an exactly marginal operator which could correspond to increasing the radius of the dual 
background. In fact, the CFT has an exactly marginal double-trace operator which is a 
product of the left and right U(l) currents, Jl(z)Jr{z). This operator breaks the Af = 2 
supersymmetry as well as some of the the extended H^-symmetries. For the N = 2 coset 
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CFT (2.7), this marginal operator changes the radius of the compact scalar in the bosonized 
formulation. In the limit of large radius, a large gap develops between the dimensions of 
typical momentum and winding operators. It would be interesting to study the effect of the 
marginal deformation on the operator dimensions for N > 2, and to see if deforming the 
CFT along this marginal direction could also create a large gap in the spectrum of operator 
dimensions. The presence of such a gap would suggest [EE] a weakly curved AdS^ dual of 
the large N CFT. 
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Appendix A: Fundamental matter 



This appendix briefly describes the high density physics of the theory in ( 1.1 ), but for the 
case where \I/ transforms as a fundamental of the SU (N) gauge group. For completeness, we 
include the case where \1/ has a flavor index which takes Nf values, and then the model has 
a U(Nf) global symmetry. It is convenient to decompose the global symmetries into a U(l) 
symmetry associated with the 'charge' density, and a SU(Nf) flavor symmetry (the latter is 
absent for Nf = 1). We proceed just as in section [Hj The high density limit is characterized 
by a Fermi wavevector of gauge-charged fermions given by 

£ = NN& (Al) 

L 71 



note that the r.h.s. has a prefactor of N, rather than the (N 2 — 1) in (1.2). The fluctuations 
near this wavevector map onto the m — // = theory, which was considered in [HTJ 188] . 
Now we bosonize the NNf complex Dirac fermions differently. Rather than considering 
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them as 2NNj Majorana fermions, we note that they can be used generate WZW currents 
of SU(N)n,, SU(Nf)N, and U(l), and these fully span the Hilbert space [87H89] . The 
SU(N) Nf currents are projected out by the gauge field, and so the low energy theory is 
made up of two decoupled sectors: a c = 1 free gapless scalar associated with the U(l), 
and a SU(Nf)^ WZW model associated with the SU(Nf) global flavor symmetry with c = 
N(Nf — l)/(N+Nf). This decoupling of the U(l) density mode is the key simplifying feature 
of the fundamental matter case, and was absent for the adjoint matter case considered in the 
body of the paper. Consequently, the U(l) sector here is similar to an ordinary Luttinger 
liquid, while the spectator SU(Nf)^ WZW model is directly linked to the additional global 
flavor symmetries of the model. 

In the notation of [HD], we can write the Hamiltonian of the U(l) sector as 



H = — I dx 

2tt 



^ (d^f + k (d x ey 



(A2) 



where K is the Luttinger parameter, 9 and <fi are scalar fields obeying the commutation 
relations 

[d x cf>(x), B{x')} = [d x 9(x), 4>{x')} = m6(x - x'), (A3) 

and the U(l) charge is 

Q = - [ dxd x <f>. (A4) 



7T _ 

The variable K is related to the exactly marginal perturbation to the U(l) theory, the analog 
of the 'radius' of the scalar in the string theory notation. The fermion fields are related to 
these scalar fields via [91] 

i> R , L = e-^vm^ kL ^ L (A5) 

where <p c R L is the primary field the SU(N)^ f WZW model transforming as a fundamental of 
SU(N), <Pr L is the primary field the SU^Nf)^ WZW model transforming as a fundamental 
of SU(Nf). The exponential factor has been chosen so that the free fermion theory without 
the SU (N) gauge field is properly bosonized at K = 1 with both WZW models conformal 
so that dim^ J = (N 2 - 1)/(2N(N + N f )) and dim[^ L ] = (Nj - l)/(2N f (N + jV»). 

To determine the Friedel exponent, we need the smallest scaling dimension operator with 
Ql = 1 and Qr = —1. Applying (A5) to the operator Tr^j^^), we can set the trace over 



the SU(N)jff WZW fields to constants [HE1EH], and the scaling dimensions of the remaining 
sectors yield 

A F = A^V^] + adiznlrf] = JL + N { "lll /y (A6) 
where K is now allowed to be not equal to unity because, in general, there will be an exactly 
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marginal interaction in the c = 1 sector. We note that the Ap — > corresponds to a 
crystalline state with broken translational symmetry; such continuous symmetry breaking 
is not possible in 1+1 dimensions, but a crystal was discussed as a mean field theory of 
baryons valid in the formal large N limit [92J. 

In the fundamental matter model, the pairing operator ipLipR can be reduced to a gauge 
singlet only for N = 2. Extending our Friedel operator argument to a gauge singlet pairing 
operator implies that we need the simplest operator with Q L = Q R — — 1, and this yields 

A P = dim[e 2ie/ v 7 ^! + 2 dimb(] = — - — + —-^ , N = 2. (A7) 

P L J [ ^ Ll NN f K N f (N + N f ) ' K J 

Finally, we note that the case N — 1 and Nf = 1 corresponds to the finite density phase 
of the massive Thirring model, which realizes the simplest Luttinger liquid. 



Appendix B: Modular invariants for coset CFT2 

Modular invariant partition functions for coset CFT 2 can be constructed once invariants 
for both the numerator and denominator CFT2's have been specified [93| 194"]. Rather than 
describing the general construction, we focus on the specific example of the (N, N; 2N) 
cosets for gauged adjoint fermions. 

The numerator CFT 2 has two copies of the SU(N)n theory, each describing (A^ 2 — 1) 
adjoint fermions. There are several options for a modular invariant, including a simple 
product of the diagonal modular invariant of each of the SU(N)n factors. However, we 
should here select the invariant that describes the situation where boundary conditions 
on the combined fermions are such that a global U(l) symmetry arises. The appropriate 
modular invariant turns out to be the diagonal modular invariant of an SO(2N 2 — 2)i 
symmetry, written as 

z SO(2N 2 -2) 1 _ |^SO(2Af 2 -2)i | 2 + |^50(2V 2 -2)i |2 + 2 |^S'0(2iV 2 -2)i |2 ^g-Q 

with T, 'v,' and 'sp' denoting the identity, vector and spinor representations of SO{2N 2 — 2)i. 
This result arises from the well known result, known as non-Abelian bosonization [95] , that 
the CFT 2 based on SO(M)i, at central charge c = M/2, describes M real fermions. 

The partition sum can be re-expressed in terms of characters of two copies of SO(N 2 — l)i, 
one for each of the groups of A^ 2 — 1 fermions 

7 SO(2N 2 ~2) 1 = I 50(Af 2 -l)i~SO(V 2 -l)i 50(V 2 -1)! ~50(V 2 -1)! |2 

+ |x v 50(7V2 - ll) xf° (7V2 - 1)l + x f°^x S v ° {N2 - lh \ 2 + 2A|xf P ° (Ar2 - 1)l X SF 0(7V2 - 1)i r , (B2) 
with A = 4(1) for odd(even).The SO(N 2 — l)x characters can each be branched into 
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characters of SU(N)n. For the NS sector characters (labels '1' and V) the results are 



SO(A^-l)! _ SU(N) N SU(N) N SU(N) N 
Al A(00...00) ' A(20...1O) ' A(01...02) 

V 50(7V 2 -1)! _ SU(N) N SU{N) N 

Av — A(io...01) A(no...011) 1°°^ 

We use Dynkin labels (lih ■ ■ ■ In-i) to tag the SU(N) representions: (00 . . . 00) is the iden- 
tity, (10 . . . 01) the adjoint, etc. (a useful reference for the group theory is [HE])- We remark 
that the SU(N) representations featuring in the NS sector satisfy the iV-ality condition 

k + 2k + ...(N - l)l N -! = mod N . (B4) 

To obtain a partition sum for the coset CFT 2 the following steps are taken. First one 
writes branching rules for the SO(2N 2 — 2)i characters into products of branching functions 
times characters of the affine algebra SU(N)2N that features in the denominator of the coset. 
Schematically 

xf (2Ar2 - 2)l =E fe AXxfW-. (B5) 

a 

The labels a take values in the integral dominant weights of SU(N) 2 n, which are written as 
Dynkin labels (M2 • • • In-i) satisfying Y^ih — We find that, for general N, the terms 



on the r.h.s. of (B5) are grouped into combinations of the form 

Aa + X\( a ) 1 A A 2( a ) + • • • l C0 J 

where A is the automorphism 

A : (kl 2 . . . l N ^) -+([2N-J2 k]hh ■ ■ ■ In-*) ■ (B7) 

i 

Again, the SU (N) representations featuring in the NS sector satisfy the iV-ality condition 



(B4) 



Writing the modular invariants for the denominator (d) and numerator (n) as 

T d 50(2V 2 -2)i_50(2V 2 -2)i 



AB 

Z» = Y,K b x! U{NhN X S b UiNhN (B8) 



ab 

the coset invariant is obtained as 



Z coset = NabN^b ■ (B9) 



ABab 
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For the (N, N; 2N) cosets, the natural choice for the denominator modular invariant 
is the SU(N)2N invariant that displays the same grouping of characters, according to the 



automophism (B7), that we observed in the branching rules (B5). Such an invariant exists 
for general N \§T\] we display explicit examples for N = 2, 3 in the main text of this paper. 

The M = (2, 2) superconformal symmetry of the (N, N; 2N) coset guarantees that the 
branching functions b a A are characters of (an extension) of the M = 2 superconformal algebra. 
For iV > 3 we find (see again the main text) that the vacuum character of this extended 
symmetry takes the form 

,A/-=2,ext i NS ip^NS , iNS , f R1 n\ 

The (N, N; 2N) coset modular invariant can be viewed as a diagonal modular invariant for 
this VT-extension of M = 2 superconformal symmetry. 



Appendix C: DLCQ Quantization of the gauged adjoint Dirac fermions in 1+1 
dimensions 

In this appendix, we will summarize the computation of the discretized light cone quan- 
tization spectrum of the 1 + 1 dimensional SU (N) gauge theory coupled to adjoint Dirac 
fermions in the large N limit. 

The first step is to write down the Lagrangian which follows the general pattern of 
[50| l52| \M\ except that the fermions are now complex. Here we follow II. B of [50]. We start 
with (1) of [54J but treat 

$ = 2 1 /* / (CI) 

as Dirac fermions. 

The light cone coordinates are defined by 

x ± = ^=(x°±x 1 ) (C2) 

so that 

A^-^oi^i). (C3) 
We will use for the Dirac matrices 
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The Lagrangian is normalized to take the form 



Tr {^(i$ - m)#]} = Tr \2iip j d+ip + 2i\ ] d^x ~ iV2m(^x + xV)) • 



(C5) 



To compare this Lagrangian with the Hamiltonian (2.1), simply note that the fermion field 
\& can be expressed in the standard mode expansion 

dkt 1 



u(k)p(k)e~ lk ^ + v(k)h(k)e l 



2tt y/2ko 

where u{k) and v(k) is the standard 2 component spinor basis satisfying 

(ft — m)u(k) = (ft + m)v(k) = 0, u(k)u(k) = —v(k)v(k) = 2m 



(C6) 



(C7) 



Then, in terms of p(k) and h(k) we recover (2.1) for the Hamiltonian. 



Returning to (C5), we gauge the free fermion theory by introducing covariant derivatives 



D* = d l Stf + i[A li ,V] . 



(Ch 



It is customary in light cone quantization to use the gauge 



A_ = 



(C9) 



so that the gauge kinetic term takes the form 



^9ym 



TrF 



2 Tr(«9_A + ) 2 



9ym 



and the Lagrangian reads 



£ = Tr | 2i^d + i> + 2i X ] d^x ~ V2im(^x + xH) + A + J + + -^—(d-A^ 
I 9ym 



with 



(CIO) 



(Cll) 



(C12) 



If we take x + as the time direction, \ an d A + are non-dynamical and can be integrated out, 
giving rise to the light-cone momentum and Hamiltonian 



P- 



dx~ Tr {2i^d-ip} 



,2,/.t_L/, - - — 



dx Tr <^ -im tp 1 ' - -g YM J + 



(C13) 
(C14) 
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Imposing canonical quantization on the fermions gives rise to relation 

^M^ - )} = ^H x ~ - %~) 0>ii6jk - ^SijSk^j (C15) 

with 

{M*-),Mz-)} = {^( a? ~)^«(*")} = ° ■ ( C16 ) 

The Dirac fermions are expanded in modes 

Hx~) = -t W E {Mn)e- mnx ' /L + Bj t (n)e—~ ' L ) (C17) 

^ neodd>0 

where we have compactified the x~ direction and imposed the anti-periodic boundary condi- 
tion on the ifj(x~) field; this typically leads to a better DLCQ computation than choosing the 
periodic boundary condition and removing the zero mode by hand. The choice of boundary 
condition should not matter in the decompactification limit K — > oo. 
The anti-commutation relation for the modes are is 

{A i:j (m),Aki(n)} = 5(m + n) (5 a S jk - -^<%4^ , (C18) 

{Sy(m),S fcI (n)} = 8(m + n) (s u 5 jk - ^SijSk^J (C19) 

where n takes odd integer values, and 

A^-n) = 4(n), Sy(-n) = £j» . (C20) 

We can now set up the light cone vacuum 

Aij(n)\0) = %(n)|0) = 0, (n > 0) . (C21) 

The states are then generated by acting by a string of "letters" A(— n) and B(—n) in a 
single trace state, i.e. 

|V) = #Tr(B(-n 1 )A(-n 2 )...B(-n k ))\0) (C22) 

where # is a symmetry factor to ensure that the norm of each state is one. 

Our next step is to write the light cone momentum and Hamiltonian operators in terms 



of the fermion oscillators. It is clear that the light cone momentum operator (C13) can be 
written in the form 

P + = J2 { T 4H^» + J-Bl^B^n)} = — , (C23) 

n>l 
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when acting on a state with fixed K. 

Instead of writing the light cone Hamiltonian P~ in terms of fermion oscillators, let us 
consider the Lorentz invariant mass operator 



2P + P- = V + T (C24) 
where V corresponds to terms associated with the term proportional to m 2 and T corre- 



sponds to the term proportional to (J + ) 2 in (C14). Then, it is easy to show that 



V = Km " E { (I) A U n )Mn) + 4(n)S i4 (n)} . (C25) 

Note that the dependence on L drops out, but there is still a dependence on K. 

Computation of T involves a somewhat tedious exercise of normal ordering the {J + ) 2 
written in terms of fermion oscillator operators. One can organize T in accordance to the 
number of oscillators destroyed and created. 

T = 2W + TUa + T 2 ^ 2 + T 3 ^i . (C26) 

In this form, one finds after some algebra, that 

= 2 jh^Y^\(— ) 2 ^(n)Mn)+(— ) 2 4(n)B,,(n)](C27) 
7r \ n — to/ J \n — m J J 

n m=l I v x » 



T^ 3 = (^^)^5(n 1 + n 2 + n 3 -m 1 )x 



1 Mk^MlMAUn^A^) 



- r4(n 3 )4(n 2 )SJ fc (^i)^K) 



- -4(n 3 )5| fc (n 2 )4(^i)^K)i (C28) 



(n 3 - mi) 
2 

7- _ m ^ kl {n2)A\ j {n 1 )A ki {m 2 )A i j{m 1 ) 

2 2 \ 

7 \o - 7 ; so A\ k (n2)Bl l (n 1 )A ij (m 2 )B j i(m 1 ) 
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(mi + m 2 ) 2 
2 



2 



(rii-mi) 2 (m l + m 2 ) 2 



B li (n 2 )AUn 1 )B ki (m 2 )A ij (m 1 ) 



3-> 1 



(ni — mi) 2 ' 
(9 2 ymNK\ 



(C29) 



V 2tt J 
2 

(rii — mi) 2 ' J 

2 2 



+ 



5(ni — mi — m 2 — m 3 ) x 
4(wi)5tt(m 3 )i < . ! (m 2 )A JJ (m 1 ) 

Ajj. (rai ) (m 3 ) B w (m 2 ) (mi ) 



(m 2 + mi) 2 (rii — m x ) 2 
2 

v2Ai,(ni)A fci (m 3 )A ij -(m 2 )S j ,(mi) 

^mi -f- m 2 j 

?t 



(rii — mi) 2 
2 



Bj k (n 1 )A ij (m 3 )B j i(m 2 )Bi k (m 1 ) 



+ 
+ 



(mi + m 2 ) 2 (rii — m x ) 5 



Bl l (n 1 )B ki (m 3 )A ij (m 2 )B j i(m 1 ) 



(mi + m 2 ) 



2 Bj i (n 1 )B jl (m 3 )Bi k (m 2 )A ki (m 1 ] 



(C30) 



At this point, a computer program must be written to generate the set of states and the 
elements of the mass operator M 2 = 2P+P-. 
As an example, for K — 5 and Q — 1, we find 



T = 



9 2 ymNK 
2tt 
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1 

4 








3 
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/ 



(C31) 



whose eigenvalues in units of gy M N /2tt are 

{0,6.25,10,10,17.5} 



(C32) 



Strictly speaking, eigenvalues of just the T without the contribution from V are unreliable 
since they correspond to taking massless fermions as the matter fields. Nonetheless, it is 
encouraging to find an exactly massless state in the spectrum which would survive the limit 
of strong gauge coupling. The massless states for these m 2 = cases continue to be present 
as the values of K are increased. 
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K = 5 


K = 7 


K = 9 


K = 11 


K = 13 


if = 15 


K = 17 


K = 19 


if = oo 


M/Q 


Q = l 


4.60 


5.08 


5.42 












6.42 


6.42 


g = 3 


3.42 


3.67 


3.85 


4.00 


4.10 








4.49 


1.50 


Q = 5 




5.51 


5.87 


6.14 


6.36 


6.54 






7.40 


1.47 


Q = 7 






7.55 


7.97 


8.30 


8.58 


8.82 




10.19 


1.46 


Q = 9 








9.57 


10.03 


10.041 


10.73 




12.83 


1.43 


Q = ll 










11.59 


12.08 


12.49 


12.85 


15.54 


1.41 



TABLE III: Mass M of the lightest fermionic bound state in the fixed Q sector for various K. The 
"— " indicates entries which are not defined. The "• • • " indicate entries which are well defined but 
have not been computed due to limits in computational resources. These are the data presented 
in figure [3} Note that M as a function of Q is minimized at Q = 3. However, M/Q as a function 
of Q appears to slowly be decreasing monotonically. 





K = 4 


K = 6 


K = 8 


K = 10 


K = oo 


M/Q 


Q = 


2.31 


2.47 


2.58 


2.67 


2.88 


oo 



TABLE IV: Mass M of the bosonic bound state in the Q = sector for K = 4, 6, 8, 10. 



The actual computation reported in section III is the computation of the spectrum of 

M 2 1 



AV + T) 

where, for definitiveness, we set the dimensionless parameter 
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The result of carrying out the calculation for K = 5,7, 9, 11, 13, 15, 17, 19, is summarized 



in table |III[ There, we tabulate the calculated value of the hadronic bound state mass M 
for fixed Q as K is increased. The K = oo is a result of linear extrapolation illustrated in 
figure |3j We observe that the lightest bound state appears to be in the Q = 3 sector, at least 
among the fermionic states. However, the quantity M/Q which determines \i crit appears to 
be slowly decreasing as Q is increased. 

Similar computations can be carried out for the bosonic bound states when the values of 
K are taken to be even. We did not perform the computation at the same level of precision 
for the bosonic bound states as we did for the fermions. The plot analogous to figure [2] is in 
figure [4] below. From figure |4j it is quite apparent that the pattern of states with smallest 
Mf Q is becoming degenerate at around \i cr n = M / Q w 1 as Q is increased. 

There is one additional feature which is notable regarding the bosonic spectrum: the state 
with smallest M is in the Q = sector. For K = 4,6, 8, 10, we find the masses presented in 
table ETJ 

A closer examination of the wavefunction indicates that the lightest state is mostly a 
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FIG. 4: Spectrum of bosonic bound states for K = 4, 6, 8, 10. This figure is the analogue of figure 
[2] for the fermionic bound states. 



mixture of the "two bit" states of the form 

^2c p TtA\p)B\K -p)\$) . (C35) 

p 

This is analogous to what was found for the adjoint Majorana model [52]. These are inter- 
esting features of our model from the point of view of its dynamics at vanishing chemical 
potential. As should be clear from the right most column in table |IV[ however, the Q = 
states do not have any impact on the physics at finite chemical potential \i. 
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